CTP #2217 
June 1993 



A Formula for the Static Potential Energy 
in Quantum Gravity. 



Giovanni Modanese^ 

Center for Theoretical Physics 
Laboratory for Nuclear Science 
Department of Physics 
Massachusetts Institute of Technology 
Cambridge, Massachusetts, 02139, U.S.A. 



Abstract 

We give a general expression for the static potential energy of the gravitational inter- 
action of two massive particles, in terms of an invariant vacuum expectation value of the 
quantized gravitational field. This formula holds for functional integral formulations of eu- 
clidean quantum gravity, regularized to avoid conformal instability. It could be regarded 
as the analogue of the Wilson loop for gauge theories and allows in principle, through nu- 
merical simulations or other approximation techniques, non perturbative evaluations of the 
potential or of the effective coupling constant. The geometrical meaning of this expression 
is quite simple, as it represents the "average proper-time delay" , respect to two neighboring 
lines, of a very long geodesic with unit timelike tangent vector. 
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1 Introduction. 



The present paper is concerned with the problem of the energy of the gravitational field. This energy has 
been under investigation since the birth of General Relativity and some issues, like the determination of the 
total energy of a field configuration, have been settled in a rigorous way in the ADM formalism or through 
Noether's theorem Q. Other points, however, like the possibility of "localizing" the gravitational energy, are still 
obscure. Before presenting our contribution, which concerns in fact the particular issue of the static potential 
energy, we shall briefly review a few general facts. 

In principle, the energy of the gravitational field is physically as important as the energy of any other field. 
In fact, one of the basic principles of relativistic field theories is that any interaction between two particles is 
not instantaneous, but it is transmitted by a field which propagates with finite velocity. Suppose that the two 
particles exchange an amount of energy as a result of their interaction. If the first particle loses the energy 
E at the time t, and the second particle receives that energy at the time t + At, it is usually assumed that in 
the interval {t, t + At) the energy is "stored" in the field that carries the interaction. 

In practice, however, the gravitational energy turns out to be much more "elusive" than other forms of 
energy, in the sense that it seems not possible to "localize" it. 

For example, the electromagnetic field is known to possess the local energy density 

r°«(x) = i[E2(:r)+B2(x)], (1) 

which is a component of the energy-impulse tensor 

^ F^^i^"; - ^ri^"'F''^Fc,0. (2) 

(See § I for our convention about the signature of the metric and others.) The tensor r'^" is locally conserved, 
that is 

d^T^'" - 0. (3) 

This assures that the spatial integral of r^" is conserved, provided E and B vanish sufficiently fast at infinity. 

Eq. (H) has suggested a generalization of the tensor t^^^ to the gravitational case, called the "Bel superenergy 
tensor" It has the form 

where JI*p°"^P is a suitably defined dual of the Riemann tensor. 

The tensor 7''"''"^ has a positivity property, is covariantly conserved and allows to define a criterium for 
the presence of a gravitational wave in a point of empty spacetime. We remind that, due to the equivalence 
principle, the curvature tensor is the true physical field in General Relativity, since we cannot set it equal to 
zero at a given point by a transformation of the coordinates. So it could be fairly expected that a local energy 
density contains this tensor. However, the tensor T^'^^'^ defined in (^) is quadratic in the curvature, and this 
prevents us from relating it to the Hamiltonian of General Relativity. 

Another problem for a local energy density resides in the fact that the ordinary conservation equation of a 
tensor (like (eq. |^)) is not generally covariant, and the introduction of the covariant derivative spoils in fact the 
conservation. All we can obtain in General Relativity is thus a conserved energy- impulse "pseudotensor" [|[ . It 
leads to an energy which is the integral of some derivatives of the metric on a surface at spatial infinity, namely 
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where it is assumed that the metric is quasi-minkowskian at infinity and is decomposed as g^v{x) = rifj_^ + hfj_^(x). 

The most natural way to obtain the pseudotensor is to apply Noether's theorem to the case of coordinates 
transformations, whose generators are the Lie derivatives The nice formalism of "improvement" of the 
Noether currents also works in this case ^J; we symmetrize the tensor obtained from translational invariance 
and we finally obtain for the energy density a quantity which amounts to a spatial divergence and reduces to 
upon spatial integration. This result is typical of the application of Noether's theorem to gauge symmetries, 
which usually produce "trivially conserved currents" of the form 

J^.=^,^, (6) 

where /^^ is an antisymmetric tensor. 

For known asympotically flat metrics like the Schwartzschild metric, (||) reproduces the total energy obtained 
from the ADM hamiltonian ^. (See also ref. 0] for a comparison of different "quasilocal energies" and 
references.) 

The purpose of this paper is to relate the potential energy of the gravitational interaction of two particles 
to an invariant quantum-field average, which might play in gravity the same role the Wilson loop plays in the 
usual gauge theories. 

In other words, this energy turns out to be related in a general fashion to the vacuum average of a simple 
and well-defined invariant functional of the field. While evaluation of this average for a weak field on a fiat 
background yields the usual Newton potential energy, non-perturbative evaluations of the same average are 
likely to give rise to modifications in the coupling constant or in the dependence of the energy on the distance 
between the particles. We shall prove, however, that the energy is always negative. 

Like in the case of the Wilson loop, our formula can be implemented quite naturally on a lattice version of 
the theory, in order to allow numerical computations. 

The outline of the paper is the following. In § mainly in order to fix our conventions, we compute the 
Newton potential starting from the graviton propagator. In § ^ we recall the connection between the ADM 
mass formula and the static gravitational potential. Through the ADM formula it is possible, in principle, to 
find the relativistic corrections to the Newton potential. In § ^ we give a formula which, treating two masses as 
external sources for a quantized gravitational field on a fiat background, allows to write the potential energy of 
their interaction. This is done using a known technique of euclidean quantum field theory |pT| . Such a formula 
would allow to compute the quantum corrections to the Newton potential. In § ^ the definition of the external 
source is generalized, avoiding use of the background metric and introducing the idea of nearby parallel lines in 
curved space. Finally, in § ^ we generalize the formula for the potential energy to the case of "strong" gravity 
and in § ^ we suggest some possible applications. 



2 Newton potential. Conventions. 

In this paper we work in (3-1-1) dimensions and follow the conventions of Weinberg Q. Our metric has signature 
(— 1, 1, 1, 1). The Einstein action is given by 

^Einst. = -7^-77 ( d'^x g{x) R{x) (7) 
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and the action of a material particle of mass m is 

5'Mat. ^ ~m dp sj-g^,^ [a;(p)]i^ (p)i'^ (p) , (8) 



where x^{p) is the trajectory of the particle and p is any parameter. The dots will always denote differentiation 
with respect to the parameter. Finally, the Einstein equations have the form 

-Rm- - \9^^uR = -SttGT^.. (9) 

Let us decompose the metric in the traditional way 

g^L^{x) = rjf^^ + hf_,^{x) (10) 

and denote the linearized Einstein equations in the harmonic gauge as 

KP^hp^^T^,. (11) 

The inverse of the kinetic operator K is the well-known Feynman-De Witt propagator 



Kf^^pAx -y) = {h^u{x)hpAy)) = — — — — ' ■ (12) 



Let us now compute the Newton potential starting from the preceding equations. The field produced by a 
generic four-momentum source Tp^ in the linearized approximation is given by 

V(^) = / d'y [K-Xli^ - y) TpAy); (13) 

when the source is a particle of mass m at rest in the origin, the only non vanishing component of Tpa- is 

Tooiy)^mS^iy), (14) 

so we have 



/ioo(x°,x) = I d^y[K-Xoi^-y)m6Hy) 



oo 



dy" / d'y 



(x — y)^ — {x*^ — y^)'^ — it 



^ y x2 - (a;0 - y0)2 - ze " |x| ' ^ ' 

This is the correct result, since in the newtonian approximation we have 

goo = -1-21/, (16) 

where V is the Newton potential. We shall encounter the integral appearing in (^5|) also in § ^. 



3 Potential Energy Versus ADM Energy. 

In classical General Relativity the total energy (mechanical -f gravitational) of a physical system is given by 
the ADM mass formula (||), which has the remarkable property of involving only the gravitational field on a 
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surface at spatial infinity. We shall briefly recall here the connection between the ADM energy and the static 
gravitational potential. 

A generic static metric g^i, can be written at spatial infinity in the form 

9^J - % + ^^M^x^xj + O (r^ j . (19) 



Performing the integral (|5|) with hij given by one sees that M2 is the ADM energy ( "total mass" ) . On 
the other hand, Mi is the mass observed by measuring the newtonian force at infinity. Substituting ( |l7| ) - ( |l9| ) 
into Einstein's equations = 0, it is easy to see that Mi = M2. This is a quite natural result In other 
words it means that, according to special relativity conceptions, the source of the newtonian field is not only the 
mass, but also the energy density. For instance, in the gravitational collapse of a star a part of the gravitational 
energy is converted into kinetic energy and eventually this energy is employed to produce heavier elements from 
hydrogen or helium. If we disregard the radiation emitted into space, the newtonian field far away from the 
star remains unchanged during the whole process and the same holds for the ADM mass, which is a conserved 
quantity. 

The gravitational potential energy can be found, by definition, assuming a static distribution of matter and 
computing the metric it generates at infinity. This has been done by Murchada and York 1^ for a spherical 
matter distribution of uniform density, using conformal transformations and a special formulation of the initial- 
value equations of General Relativity. For a sphere of (small) density p and unit radius they found the right 
newtonian gravitational binding energy, namely the ADM mass is given in this case by (reintroducing the radius 
R and the velocity of light c) 

Mtot = ^ttR'p - li^7r2p2Gi?5 + Qip""). (20) 
6 15 

Remembering that AfxoT also represents the effective source of the newtonian field, we see that the second 
term in (^0|) gives rise to a deviation from the famous law which states the independence of the potential on 
the radius of the source. Nevertheless, this effect is usually unobservable, due to the very small factor c^^. 

It is also possible to find the following corrections to (pO|), proportional to p"^, p**, ... They denote the existence 
of general-relativistic corrections to the potential energy mim2G /r. For instance, the term proportional to 
would contribute to Mtot a term of the form 

AAf cx ^p^G^i?^ + 0(/). (21) 

In the case of a source constituted by two pointlike bodies of masses mi and m2, kept at rest at a fixed 
distance r, the method of solution mentioned above is not applicable. (An approximate solution is still possible 
1^.) From the preceding discussion we may infer that the ADM mass is given in this case by 

1 Gmim2 / 1 \ 
Mtot ^ mi +1712-— — + — (22) 



We are not able, however, to deduce the relativistic corrections to the two-body potential from (|2l|), because 
the corresponding potential does not admit a continuum limit. For instance, if we try to integrate a potential of 
the form G'^rrC'i'^rrC'^'^ jr'' to obtain the term proportional to p^, we find that the binding energy of the sphere 
depends on the way it has actually been put together. 
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4 Quantum formula for the potential energy on a flat 
background. 

The same result we found in the preceding Section using the classical equations of motion can be obtained 
in a completely different way. It is known that the ground state energy of a system described by an action 
So[<j)] = J d^x L{(p{x)) in the presence of external sources J{x) can in euclidean quantum field theory be expressed 
as 

^_ J exp {- / d^x Ljcpjx)) + / d^x cj>{x)J{x)) 
^-T^oo r'°S /dMexp{-/d%L(</)(a;))} ' ^'^^ 

where, outside the interval {—\T, the source has been switched off. 

This formula has been proved exactly in perturbation theory [ pT| for the case of a linear local coupling 
between the field and the external source, but it can be generalized if we assume that in any case the vacuum- 
to-vacuum transition amplitude is given by 

_ /rfHexp{-5'o[0] + 5'intcr. [0, J]) , . 

^ ' /dMexp{-5oM} ■ ^ ^ 

In fact, inserting a complete set of energy eigenstates we can write 

(0+10-). = (0|e-«^|0) = Y.^Q\e-"^\n){n\Q) ^ J] |(0|n)p e'^^^ . (25) 

n n 

The smallest energy eigenvalue En corresponds to the ground state, and in the limit T — > cx) it dominates the 
sum. So taking the logarithm and dividing by (— T) we obtain that energy. This is a well-known technique in 
QCD (see for instance Jl2t). 

In the case of a weak gravitational field quantized on a flat background, we may consider the source consti- 
tuted by two masses mi, m2, placed at rest near the origin at a distance L each from the other (see eq. (p7|)). 
The action of this system is 



S — S'Einst. + -SMat.,! + <S'Mat.,2 



^ J d'^x^M^)R{x) - mi J'^ dti ^J-g^^[x{t^)]x^^{t^)x'^{t^) 



-7712 



dt2 ^J-9^.Ay{t2)]r{t2)rit2), (26) 



where the trajectories a;^(ii) and y^{t2) of the particles with respect to the background are simply given by 

= (h. 0, o) ; yni2) - fi2, ^, 0, o) . (27) 



2' ' V " V 2 

So we have, denoting by S the euclidean action, 

£ = lim X 

T-.00 T 

J d[h] exp I -S'Einst. - J dti yj\ - /ioo[a;(ii)] - ma / dt2 \J\ - /ioo[y(^2)] \ 

log ^ 7—. ^ 

J d[h] exp I -S'Einst. j 

(28) 
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Since it is known that the euchdean Einstein action is not bounded from below, due to the existence of 
"conformal singularities" (see [p^ ), we actually mean by S'Einst. a regularized version of (|^). Such a regularization 
can be achieved adding i?^-terms to the original Einstein lagrangian jl^ or through a more recent technique, 
called "stochastic regularization" |13 . We shall return to this point in § |^. 



Returning to (28), by standard perturbation techniques it is straightforward to see that for weak fields and 



to lowest order in G it reduces to 

fyi-yfyi2 f ^ f ^ I 

IH — t'^^^ T {^oo[x{ti)]hoo[yit2)]) > 

X. n 

~ TOi + TO2 + lini - — — / dti / dt2 ^oooo('^i ~ "^2, L, 0, 0) 

= mi +1712 -T^- (29) 

The next term in the perturbative series is the first quantum correction, proportional to fiG^ . We shall compute 
it exactly in a forthcoming paper. 

Eq. (|2^), like the corresponding ones in QED or QCD, has the physically appealing feature of showing how 
the force between the sources ultimately arises from the exchange of massless bosons. However, let us make a 
closer comparison with electrodynamics. In that case the analogue of the functional integral which appears in 



the logarithm of (|28|) has the form |10 




explg dhAo[x{h)] -g ' dt2Aa[y{t2)] \ ) . (30) 



(The two charges have been chosen to be opposite: qi — g, q2 = ^.9-) Reversing the direction of integration in 
the second integral and closing the contour at infinity, one is able to show that the quantity (|30| ) coincides with 
the Wilson loop of a single charge g, thus giving a gauge invariant expression for the potential energy. 

In gravity this is not possible: we may imagine that an expression like (|3^) could be obtained in the first-order 
formalism (with Aq replaced by the tetrad Eq), but the masses necessarily have the same sign. 

Luckily, our formula for the energy is invariant as it stands, as we shall see better in the next Sections, where 
we generalize it to the case where no background metric is fixed. 



5 Parallel lines in curved space. 



In order to generalize eq. (28) beyond the case of weak fluctuations of the gravitational fleld around a fixed flat 



background, we need a definition of the source that does not depend on such a background. 
We assume that a functional integral for regularized euclidean gravity exists, denoted by 

z ^ J d[g] exp I -S'Einst. [ff]} , (31) 

and we require that all the field configurations in this functional integral are asymptotically fiat. 

Let us suppose that a field configuration is given. We consider a geodesic line of length T , which starts at 
an arbitrary point in the "past" asymptotically flat region with unit timelike velocity. 
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To fix the ideas, this curve could be written in its first part as 

e{T)= (^-^T + T, 0,0, 0<T<f, (32) 

where r is the proper time measured along the curve and we have chosen the spatial coordinates of the starting 
point to be equal to (0, 0, 0) (this is an irrelevant arbitrariness, since at the end we shall integrate over all the 
configurations of the field). As usual, T denotes a very long time interval. After a time ~ f the curve enters 
the region of spacetime where the gravitational field is non vanishing. It continues as a geodesic, which means 
that ^^(r) satisfies the equation 

r^.KM]e'(r)r(r)+|''(r)=0, (33) 

where is the ChristofFel symbol of the metric. The curve terminates at r = ^T, again in the flat region. 

Let us then take in the initial point ^'^(0) a unit vector q^{0), orthogonal to ^'^(0) (for instance, in our 
preceding example, q^(0) = (0, 1, 0, 0)), and define a vector g'^(T) along the curve ^^(r) by parallel transport 
of g^(0). We remind that C^(t), being the tangent vector of a geodesic, is parallel transported along the geodesic 
itself, and that the parallel transport preserves the norms and the scalar products. Then the following relations 
hold along the curve 

i''{r)t{r)g,Mr)] = -1; (34) 
9''(r)g"^(r)g^,[^(r)] = 1; (35) 
eW{T)g,Aar)] = 0. (36) 

Next we consider two masses mi, m2, and a length L which we may regard as infinitesimal, compared to 
the scale T. We assume that the two masses follow the trajectories x^(t) and y^ir), respectively, given by 

xf^ir) = e{r)-Liq'^{T); (37) 
y^(r) = C{r) + L2q'^iT), (38) 

where Li and L2 are two positive lengths such that 

Li + L2 — L and — miLi + m2£2 = 0. (39) 

The physical meaning of the preceding geometrical construction is apparent: it represents an observer which 
falls freely in the center of mass of the system composed by mi and m2, while holding the two masses at rest 
at a distance L each from the other. This is a generalization of the source introduced in § || that is naturally 
dictated by the equivalence principle. 

We notice that if the two masses were allowed to fall freely in the field, they would not keep at a constant 
distance from each other. In fact, as it is well known from the so-called geodesic deviation equation, the distance 
between two neighboring geodesies varies according to the sign of the curvature in the region they are traversing. 

We can reparameterize the two curves x^{t) and y^{T) introducing their proper times ri and T2, respectively. 
The ratio between the proper time ri and the proper time r is given by the equation 

dTi = ^ -g^^[x{Tl)]x^'{Tl)x^{Tl)dT, (40) 

where ti = ti(t); using (|37|), (|38|), we have 

(^^)' = l + Liq'^{T)d^g,MrM'[r)tir) + 

+Li g^Mr)] {i''{T)q'^{r) + r(r)g'^(r)} + 0{LI). (41) 
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An analogous relation holds for T2. We agree to adjust the fmiction ti{t) in such a way that ti(0) = 0. Then 
we shall denote Ti(-iT) = -^T{ and ti{^T) = ^T{' . For flat geometries we have T{ ^ T{' = T. Analogous 
relations hold for T2. 

We notice that eq. takes a much simpler form in the coordinate system where ^^(r) defines one axis 
of a normal coordinates system (if they can be globally defined). In that case, we have ^^(t) — (1, 0, 0, 0), 
g^(r) and ( pl| ) reduces to 

= l + iiai5oo(T, 0, 0, 0) + O(i?), (42) 

which obviously reminds us of the equations of § ^ However, we are interested only in coordinates-independent 
quantities in the following. 



6 General formula for the potential energy. 

According to the discussion of the preceding Section, in the absence of a background eq. (^) must be rewritten 



as 



^lim log /exp \ -mi j ^ dri \J -g^u [x{ti )]xf^ (n )i'' (ti ) 



-TO2 



\ dT2 ^ -9^,,[y{T2)]y^{T2)r{r2)\) (43) 



r™<L ^exp |-imi {T[ + T") - {T^ + 



S'Einst. 

)\)^ , (44) 

5Einst. 



where, for brevity, we have denoted by brackets the functional average weighted by the exponential of the 
(regularized) Einstein action. 

Now we exploit the property, characteristic of timelike geodesies in a Lorentzian manifold, of having maximal 
length with respect to neighboring lines. This means that 

imi(r; + Tn - T{l-<5i(Li,[g])}; (45) 

im2(r^ + r^') - T{l-52{L2A9m, (46) 

where 81^2 are small positive adimensional functionals of the geometry 5, which also depend on L1.2 and thus, 
- through eq. ( ^9|) - on L and mi, m2. So we have 

£ = hm -ilog(e-^('"i+'"^)exp{T(mi(^i(Li,[5]) + m,2<52(L2,b]))}). 

~ mi + m2 - (mi(5i (Li, [g]) + m2(52 (^2, [g]))5^.^^^ ■ (47) 

The quantity in the bracket represents the energy of the gravitational interaction. More exactly, it represents, 
apart from mi, TO2, the ground state energy of the system constituted by the source described in § |[ coupled 
to a quantum gravitational field. By construction, £ is invariant with respect to coordinate transformations. 

We see from eq. (^ ) that the interaction energy is always negative. Apart from this, (|47| ) alone does not give 
us any precise indication on the dependence of this energy on the masses, on L and on G (which is contained 
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in the action). This dependence is a nontrivial result of the dynamics, as we have seen aheady in the simple 
perturbative example of § ^. We can apply to the functional integral appearing in ( ^7|) other approximation 
techniques, like the semiclassical approximation, or discretize it and use numerical methods. 

In any case, the geometrical meaning of ( ^7[ ) is quite simple. In practice, this formula implies the following: 

(1) trace a geodesic with unit timelike tangent vector through one (asympotically flat) field configuration; 

(2) measure the total "delay" of this geodesic with respect to two neighboring lines; (3) average on many 
configurations. 

7 Possible applications and concluding remarks. 

The functional integral representation for the static gravitational potential given in ( ^ ) can serve as a basis for 
various approximations (weak-field, semiclassical) and for numerical simulations. The latter require a discretized 
version of the theory in which the geodesic lengths or the metric represent the fundamental variables. The 
most suitable technique under this respect is probably the "quantum Regge calculus" of Hamber Q ; another 
candidate, which might have the advantage of simpler algorithms, is the "stochastic stabilized gravity" of 
Greensite [ p3[ . As we pointed out in the last Section, the geometrical meaning of our formula is quite simple. 
This should allow us to write a dedicated algorithm for its evaluation, in order to improve the efficiency of the 
method. 

The most simple quantity to be "measured" in this way in a lattice version of gravity is the effective coupling 
constant G. In order to compute it through £, we just need to set mi = TO2, Li — L2 and assign some fixed 
values to them, so that the algorithm can be even more simplified. 

There have been some suggestions (see for ex. ||l^) that the effective coupling constant in gravity is scale 
dependent and decreases at small distances. This "antiscreening" of the gravitational interaction seems to be 
quite natural, since the longer is the cloud of virtual particles, the stronger is the gravitational force. In a 
forthcoming paper we shall employ our formula to compute the correction of order h to the potential energy, 
and possibly to exibit such an effect. 

In conclusion, we remind that this paper is part of a program which aims to study physical observables 
in four-dimensional quantum gravity. Other observables considered were the vacuum correlations at goedesic 
distance [|l6| and the loops of the Christoffel connection . 
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